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Abstract 



We use a boost-invariant one-dimensional (cylindrically symmetric) fluid dy- 
[ namics code to calculate e + e~ production from p° and u decay in the central 

o ' 

rapidity region of a central S+Au collision at yfs = 20 GeV/nucleon. We 

O ' 

\ use equations of state with a first-order phase transition between a massless 

pion gas and quark gluon plasma, with transition temperatures in the range 



150 — 200 MeV. The production cross section at the p mass loosely constrains 
the transition and freeze-out temperatures, and we find that the mx spectrum 



is a good thermometer for sufficiently high T c . 
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One of the most important problems in nuclear physics today is the study of the equation 
of state (EOS) of high-temperature and high-density hadronic matter; the former is necessary 
to understand the first few microseconds of the big bang, and the latter for elucidating the 
physics of neutron stars and similar dense astrophysical objects. In order to investigate 
the nature of the transition from normal (T = 0) hadronic matter to deconfined quark 
gluon plasma (QGP), or whatever the high temperature phase may be, experiments are 
being performed at Brookhaven's Alternating Gradient Synchrotron (AGS) and CERN's 
Superconducting Proton Synchrotron (SPS), and future experiments are being planned for 
Brookhaven's Relativistic Heavy Ion Collider (RHIC) and CERN's Large Hadronic Collider 
(LHC). In these experiments, the only particles expected to emerge from the early part 
of the collision (and thus the only probes of the high-temperature EOS) are photons and 
leptons, which therefore may provide the most important information obtained from these 
collisions. 

Predicted e + e~ transverse mass, %, spectra from resonance decays in ultra-relativistic 
nuclear collisions are sensitive to the hadronic transition temperature, T c , when transverse 
expansion is neglected Thus, measurements of these spectra may place strong con- 

straints on T c , as well as other parameters of the hadronic EOS ||. However, the photon 
transverse momentum spectrum was also predicted to be sensitive to T c Q, but the sensi- 
tivity disappeared when transverse expansion was included ||. The next natural step is to 
study e + e~ transverse mass spectra with a model including transverse expansion. This is 
the subject of the present paper. Since transverse expansion may increase the contribution 
of mesons decaying after freeze-out, a careful treatment of freeze-out is necessary. We find 
that the m-r spectra are in fact dominated by the freeze-out contribution for sufficiently 
high ttiti but can still be fitted with apparent temperatures correlated to T c . This may pro- 
vide a method to determine T c experimentally, even though the underlying physics is more 
complicated than originally envisioned in the treatment of lepton pairs without transverse 
expansion [|TJ. We show that the correlation between T c and the rrij- spectrum of dileptons 
from p° and uo decay persists in the presence of transverse expansion, and thus that a deter- 
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mination of the transverse mass spectrum of dileptons in the p° — 00 peak would yield useful 
constraints on T c . 

Here we adopt the position that in strongly-interacting matter at sufficiently high tem- 
perature and/or density a transition takes place to a high-energy phase. In the present 
paper, we are interested in the possibility of extracting the transition temperature assum- 
ing a first-order phase transition. We assume a boost-invariant longitudinal expansion as 
discussed by Bjorken ||, coupled to a cylindrically symmetric transverse expansion. For 
the initial evolution, we use a simple model for production of longitudinally free-streaming 
hot matter [|7| . After the hot matter has thermally equilibrated, we use thermal equilibrium 
fluid-dynamical evolution, but consider deviations from chemical equilibrium in the high- 
temperature phase by allowing the quark and antiquark densities to be a (fixed) fraction, x, 
of their equilibrium values. When particle mean free paths become comparable to the radius 
of the cylinder of hot matter, which we take to occur at some freeze-out temperature, Tf a , 
we assume that the particles crossing the freeze-out surface T = Tf Q stream freely until they 
either decay or reach the detectors. We then calculate the e + e~ production rate from p° and 
uj decays; these are not separated, as the masses are almost degenerate and the experimental 
resolution is typically not good enough to resolve the individual peaks. We investigate the 
sensitivity to different assumptions about the initial temperature, freeze-out temperature, 
and quark fraction, and compare production rates to the preliminary NA45 data |§. We 
use standard high-energy conventions, c = % = ks = 1. 

Here we describe the initial conditions and assumptions about the EOS; the details of the 
fluid-dynamical calculation can be found elsewhere |9]||. We consider a central collision of 
two large nuclei at SPS energy (-y/i = 20 GeV/nucleon). For such high collision energies we 
expect approximate longitudinal boost invariance ||, so the behavior of the produced matter 
at different rapidities is the same in the longitudinally comoving frame for fixed proper time 
t = \Jt 2 — z 2 , where z is the distance along the beam axis. At r = the colliding nuclei reach 
the point of maximum overlap and are assumed to form a longitudinally expanding pancake. 
The hot matter has thermalized at r = tq (~ 0.2 fm/c) |To|Jl~T|| , when a cylindrically- 

3 



symmetrical transverse expansion begins, coupled to the longitudinal expansion. 

From r = until the transverse expansion starts at r = r , we assume a boost-invariant 
cylinder of radius R < (the radius of the smaller nucleus), filled uniformly with QGP at 
temperature T = Tq. This is approximately compatible with the initial entropy density for 
short times [[/J. We determine T by assuming entropy conservation for r > r , hence 

3.G dN^/dy 

7T K%. Tq 

where s is the entropy density, with total (charged plus neutral) multiplicity density dN^/dy. 
The equations of state (EOS's) that we use here are of the form 
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where g q is the number of massless degrees of freedom in the deconfined phase. We treat 
only the case of zero baryon density, so the entropy density is s = (e + P)/T independent of 
the phase of the matter. Below T c , the EOS is that of a massless pion gas. Because recent 
calculations have predicted that the quarks may reach only a fraction of their equilibrium 
number by the beginning of transverse expansion |10,ll|l, we take g q = 16 + 21a;, where x is a 



parameter that we vary between and 1 to simulate the effect of reducing the quark density 
in the QGP below the equilibrium value (x=l is equilibrium for two flavors of massless 
quarks). The vacuum energy density in the deconfined phase, B, is related to the transition 
temperature, T c , by requiring equal pressures in the deconfined and hadronic phases at 
T = T c so that B = it 2 (g q — 3)T 4 /90. The value of B we use away from equilibrium is 
calculated by assuming x — 1 and supplying an equilibrium transition temperature T' c . It is 
this T' c that is actually used on the graphs and is useful for comparison purposes. 

We calculate the central rapidity region distribution from both interacting and frozen- 
out hot matter. The contribution from the interacting hot matter is 
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Here ttit = \Jm 2 + pf,, i? is the e + e~ production rate per unit four-volume in the fluid 
frame, v is the transverse velocity of the fluid in the cell characterized by proper time r, 
space-time rapidity r\ and radial position r (measured in the frame moving with transverse 
velocity zero and longitudinal velocity tanhr/ in the lab), and 7 = (1 — v 2 )^ 1 / 2 . The e + e~ 
production rate from p° and u mesons in thermally and chemically equilibrated hadron gas 



is 



dR m(g p oYpO_ 
d 3 p 
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where g p o — — 3 are the degeneracies of the p° and uj mesons, and the partial widths are 

- e - = 0.60 keV. E is the energy measured in the fluid frame, 



6.77 keV and 17 



and we use the average mass of the p and uj, m = 0.775 GeV. No mesons exist in the QGP, 
and hence the resonant e + e~ production rate is zero for that part of the fluid. 

Mesons which pass through the freeze-out surface T = Tf Q are no longer subject to 
fluid-dynamic flow, as their mean free paths are so long that they do not interact with 
the surrounding matter. As these frozen-out mesons free-stream toward the detectors they 
decay to e + e~ pairs as before, but the number of mesons decreases exponentially in time, 
with time constant equal to the free space meson decay time. The contribution from the 
decay of frozen-out p° and uj mesons is thus 



mTdrriTdy 



(It t 2 tt i'( t) I ' d(cos9) J** d(f) (5) 
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where 
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p sin 9 cos cf) + vy/p' 2 + m' 2 



(6) 



\Jp 2 + m 2 + sin cos 



is the component of the meson velocity perpendicular to the freeze-out surface, and T p = 151 
MeV and = 8.43 MeV are the total widths for the p and u. dr/dr is the radial velocity 
of the freeze-out surface and r(r) is its location. The G function ensures that only mesons 
going through the freeze-out surface in the outward direction are counted. 

Finally, we extract the fit temperature, Tf it , by performing a least-squares fit to our 
calculated spectrum (equilibrium plus freeze-out) in the window 1.155 < < 1.755 GeV, 
which is dominated by freeze-out. The part of the spectrum that is dominated by the 
equilibrium contribution is much smaller, typically 0.8 < < 0.9 GeV. Unfortunately this 
is of the order of a single experimental bin, so we choose instead to focus on the freeze-out 
contribution, in anticipation of fitting experimental data with both equilibrium and freeze- 
out contributions. Furthermore, a measurement of an tjit spectrum from the p — uj peak 
would look at a finite range of m, so we filter the simulated data through an invariant mass 
bin of size 0.475 < m < 1.075 GeV, assuming a Breit-Wigner distribution for each species. 
For our fitting function, we take eq. (13) with v = (i.e. ignoring transverse expansion), 
resulting in the formula ffl 



Our "standard" parameter set is r = 0.2 fm/c, x=l, Tf Q = 120 MeV, and 150 <T C < 200 
MeV; if different values of any of the parameters are given, all the others are held at the 
standard values. We use dN^/dy = 188 for central S+Au collisions at SPS energy (the NA45 
experimenters estimate dN c h/di] = 125 for their S+Au collision sample, and we assume 
isospin symmetry). Our standard value for the equilibration time, r = 0.2 fm/c, implies 
an initial temperature T = 327 MeV. Fig. 1 shows the equilibrium, freeze-out and total 
contributions for the standard set with T c =170 MeV. 
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In the absence of transverse expansion, Tfu is a monotonically increasing function, and 
comparable to T c , so one can infer T c given Tf it from the measured e + e~ spectrum [J]. In 
Fig. 2, we vary T c from 150 to 200 MeV and calculate the resulting values of Tfu while 
including transverse expansion. We find that Tfu is higher than in the case of longitudinal- 
only expansion, but still monotonic (except where Tf Q approaches T c ). Tf it is 10-20 MeV 
lower for tq=1 fm (To=192 MeV), which we attribute to the fact that transverse expansion 
develops relatively little in this case. 

Note that freeze-out is treated dynamically in the present calculation in contrast to an 
earlier investigation Jl| where transverse expansion was neglected and thus freeze-out was 
approximated as occurring at a fixed proper time. We consider a moving freeze-out surface 
here, and the yield from mesons crossing this surface is integrated over the history of the 
collision. This contribution is enhanced by transverse expansion; the larger the temperature 
difference T c — Tf D , the stronger the enhancement. We find that increasing Tf a to 140 MeV 
from its standard value decreases Tf it . For larger Tf Q the hadronic matter near the freeze-out 
surface expands outward at a lower velocity, so that there is less contribution to from 
the fluid motion. Thus Tfu, which depends strongly on the fluid motion, is lowered. 

The curve for x = 0.7 shows that Tf it is about 15 MeV higher than for x = 1. The increase 
occurs because the temperature of the mixed phase, which makes a large contribution to 
the e + e~ spectrum, is raised when x is lowered and B is held constant. 

Our overall production at the p° — u peak is up to a factor 3 below the value mea- 
sured by the NA45 experiment. We take the value measured in the experimental bin 
0.7 < m < 0.8 GeV (which contains the p° — u peak), assuming that the signal is all 
from meson decay, obtaining d 2 / dmdr] = 1.7 x 10~ 3 GeV" 1 . We calculate the same 
quantity for our simulated data from the standard calculation with T c =150 MeV, includ- 
ing the experimental acceptance by counting only e + e~ pairs whose members both have 
2.1 < i] < 2.65, and integrating over a Breit-Wigner distribution for each species, giving 
us d 2 iVg+gl 1 ^ / drndr] = 5.9 x 10~ 4 GeV -1 for our standard parameter set. Increasing T c to 
200 MeV makes d 2 N { £+" ] / drndr] = 1.2 x 10~ 3 GeV" 1 , still below the preliminary NA45 
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data. Additional pair-producing processes not included here will increase the cross-section. 
Also, one likely effect of a departure from chemical equilibrium in the hadronic phase is an 
excess of p mesons, increasing the strength of the peak and thus accounting for some of the 
missing signal. The overall production rate is sensitive to Tf a : increasing T/ G from 120 to 
140 MeV while keeping T c = 150 MeV gives d 2 N { f e t w) / dmdr] = 8.5 x 10~ 4 GeV" 1 (T fo = 140 
MeV with T c = 200 MeV gives d 2 / dmdr] = 1.5 x 10" 3 GeV" 1 ); on the other hand, 
sensitivity to increasing tq or decreasing x is minimal. Omitting transverse expansion gives 
a value that is a factor of 2-3 too high. Since the estimated error in the NA45 data is of 
order 50%, we would not rule out any of the calculations except that with no transverse 
expansion. 



Partial restoration of chiral symmetry in hot and dense matter [|12j is expected to modify 
the masses and lifetimes (widths) of vector mesons. Changes in the p-meson properties 
may indicate chiral restoration through modified e + e~ production. However, there is no 



consensus at present on the behavior of the p [IT]]. Arguments based on QCD sum rules (13 



and effective Lagrangians [T41] point to decreasing meson masses with increasing temperature 



and density, while vector-meson dominance studies |T5[ and consistency arguments JTB] have 
been used to support the opposing view. For the time being, we therefore use the free mass 
and width of the resonances as a first approximation. The study of the effects of in-medium 
modifications on e + e~ production from thermal resonances as a probe for chiral restoration 
is left for future work. 

It is tempting to carry out a calculation without a transition to a high-energy phase 
in the present model. However, this would not be sensible with the massless pion gas 
approximation. The complete resonance spectrum would have to be taken into account and 
such an investigation is beyond the scope of the present paper. 

We have shown that transverse expansion does not destroy the correlation suggested in 
Ref. If T c is close to the high end of the interval considered here (T c ~ 200 MeV) and 
Tf is not too low, we come close to reproducing the preliminary NA45 data, so that the 
shape of the e + e~ ttlt spectrum, as parametrized by Tfu, should make a good thermometer 
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to measure T c in this region. If, however, T c ~ 150 MeV or Tf < 140 MeV, the contribution 
from other (potentially long-lived) resonances appears to be significant, and therefore the 
usefulness of the p meson as a thermometer in that region is questionable. The preliminary 
NA45 data appear to rule out models in which transverse expansion plays no role, but 
otherwise place no significant constraints on the collision dynamics. 
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FIGURE CAPTIONS 

Fig. 1: A central S+Au collision at SPS energy, using our standard set of parameters and 
T c =170 MeV. 

Fig. 2: Tf it for various parameter sets as a function of T c . 
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Neumann fig 1 
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Neumann fig 2 
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